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ABSTRACT 

We derive the rates of capture N of main sequence turn off stars by the central massive 
black hole in a sample of galaxies from Magorrian etal. 1998. The disruption rates are 
smaller than previously believed with N ~ 1CP 4 — 10~ 7 per galaxy. A correlation be- 
tween N and black hole mass M is exploited to estimate the rate of tidal disruptions 
in the local universe. Assuming that all or most galaxies have massive black holes 
in their nuclei, this rate should be dominated by sub-L* galaxies. The rate of tidal 
disruptions could be high enough to be detected in supernova (or similar) monitor- 
ing campaigns — we estimate the rate of tidal disruptions to be 0.01 — 0.1 times the 
supernova rate. We have also estimated the rates of disruption of red giants, which 
may be significant (N ii 10~ 4 y _1 per galaxy) for M <; 10 8 Mq, but are likely to be 
harder to observe — only of order 10~ 4 times the supernova rate in the local universe. 
In calculating capture rates, we advise caution when applying scaling formulae by 
other authors, which are not applicable in the physical regime spanned by the galaxies 
considered here. 
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1 INTRODUCTION 



A number of early type galaxies and spiral bulges are 
now thought to contain massive black holes in their nuclei 
(Magorrian etal. 1998). Direct evidence is also now available 
which supports the idea that active galaxies are powered by 
massive black holes (Reynolds & Fabian 1997, Fabian etal. 
1998). It can be argued that many if not all galaxies could be 
expected to have undergone an active phase and to possess a 
massive black hole (Haehnelt & Rees 1993). Magorrian etal. 
measure black hole masses in the range M = 10 6 — 10 9 Mq, 
and report that M is correlated with the mass M of the 
underlying hot stellar system, with M/M — x ~ 0.006. 

A main sequence star of roughly solar type can be 
tidally disrupted by a black hole with mass M ^ 2x 10 8 Mq 
(Hills 1975). Larger black holes swallow main sequence stars 
whole, but red giants are susceptible to tidal disruption be- 
cause of their lower density. Frank & Rees (1976) (FR), 
Young, Shields & Wheeler (1977) considered a system com- 
posed of a black hole and an isothermal sphere of stars and 
derived analytic expressions for the rate of capture of stars. 
Cohn & Kulsrud (1978) (CK) were able to calibrate FR 
using a more detailed numerical calculation. Estimates of 
capture rates in the literature are often based on equation 
(66) of CK (e.g. Cannizzo, Lee & Goodman 1990, Kormendy 
& Richstone 1995). We argue that such estimates are often 
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wrong because the physical conditions in galactic nuclei are 
different from those anticipated by CK (who were in any 
case more concerned with globular clusters). 

The next section summarises briefly what is known 
about the observational signatures of tidal disruption. Sec- 
tion ^| establishes notation by reviewing the analytic theory 
of FR and sets out our method for calculating capture rates 
in the Nuker galaxies. We calibrate our calculation against 
the result of CK. In Section ^ we describe the data and list 
the quantities we derive from theim including capture rates. 
We discuss our results in Section 



2 OBSERVABILITY 

The observable consequences of such an event have been dis- 
cussed by Rees (1988), Evans & Kochanek (1989), Canizzo, 
Lee & Goodman (1990) and Ulmer (1997). The general ex- 
pectation for the disruption of a main sequence star is a 
flare with bolometric luminosity near the Eddington limit. 
In comparison, the luminosity from the central r < 1 arcsec 
of all the Nuker galaxies is small compared to the Edding- 
ton luminosity of the black hole (^ 10~ 4 Z/Edd, Table [l]). The 
bolometric luminosity should be close to the Eddington limit 
because of the predicted mass accretion rates are above the 
Eddington accretion rate for black holes with mass up to 
about 3 x 10 7 Mq . Above that the predicted mass accretion 
rate divided by the Eddington mass accretion rate falls as 
M 3//2 (Ulmer 1997). There is much uncertainty as to how 
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optically bright a tidal disruption event will be, in optical, 
but the minimum luminosity expected for typical disrup- 
tions of main sequence stars is 10~ 3 — 10 _2 LEdd for the U 
and V bands (Ulmer 1997). Typical durations of the bright- 
est phase will be 0.1 — 2 years depending on the black hole 
mass. Main sequence tidal disruption events should there- 
fore be relatively easy to detect as a byproduct of supernova 
searches (Section |4.l| ) or other searches which observe the 
same galaxies many times over (e.g. Southern strip of the 
SDSS). 

Red giant disruption events are generally much longer 
lasting, and therefore fainter, than main sequence events. 
Red giants typically are disrupted by growing onto the loss 
cone (Section fc.2| ), so at disruption, they have pericenter 
equal to the tidal disruption radius. Consequently, an aver- 
age red giant (with R ~ 100-Rq) disruption has a approxi- 

1/2 

mate duration of 100M 6 years and a produces a peak mass 
accretion rate no larger than ~ 10 -3 M 6 ~ 1/ ' 2 Mq yr _1 (see 
Ulmer 1997, eqs 3,5). The Eddington accretion rate, which 
would supply enough mass to provide the Eddington lumi- 
nosity, is 3 x 10 -2 Mee^ 1 , where 0.1 x e.i is the efficiency of 
rest mass to light conversion in the accretion process. This 
means that red giant disruption could keep a small, 10 6 Mq 
black hole illuminated at (l/30th) its Eddington rate for up 
to a thousand years. 

Just how bright a red giant disruption would appear 
in optical bands is difficult to determine, but a simple es- 
timate is as follows. If most of the energy is released from 
within a few Schwarzschild radii for both the main sequence 
and red giant disruptions, and the bolometric luminosity of 
the red giant disruptions is 1/30 that of the main sequence 
disruptions for a 10 6 Mq black hole, then 



1 

30 



Lrg 
Lms 



ff,RG 



-"-off, MS 



(1) 



Therefore, the ratio of temperatures is (1/30) 1//4 ~ 0.4. Be- 
cause the optical bands lie in the Rayleigh- Jeans tail of the 
spectrum where the luminosity scales linearly with T, we 
would at first guess expect that red giant disruptions would 
be nearly (i.e. ~ 0.4 times) as bright as main sequence 
events. Around more massive black holes with M ~ 10 s Mq, 
red giants may be significantly less luminous if the accre- 
tion becomes advection dominated as appears to be the case 
for many extra-galactic black holes with accretion rates far 
below the Eddington accretion rate (e.g. NGC 4258 (La- 
sot a etal. 1996), M87 (Reynolds etal. 1996)). Around very 
massive black holes with M ~ 1O 1O M0, the timescale for 
return of the disruption material to the black hole becomes 
comparable with the interval between successive tidal dis- 
ruptions (see Figure ^). Red giant disruptions would then 
provide a small, nearly constant flow of material onto the 
black hole of ~ lO _3 M0yr _1 corresponding to a luminosity 
of ~ 10- 5 e.iL Edd 



3 THEORY OF CAPTURE RATES 

Consider a spherical cluster of stars of mass m„ with density 
p(r), and isotropic velocity dispersion a(r). Typically p(r) is 
in the form of two power laws, separated by a break radius 
Th. The velocity dispersion is then given by Jeans' equations: 



dpa 2 GpM(r) 
dr r 2 



= 



(2) 



where M(r) is the mass enclosed within radius r, including 
the black hole, and G is Newton's constant. The black hole 
exerts an influence on the stars inside a sphere of influence 
with radius r a given by 

GM 



er(r)'=0, 



(3) 



where M is the black hole mass. Within the sphere of influ- 
ence (r < r a ) 



GM 
(l+p)r' 



(4) 



where the logarithmic slope of the density is — p (i.e. p ~ 
r~ p , with p possibly dependent on radius). An important 
derived quantity is the two-body relaxation timescale, given 
by 



tr 



AG 2 pm, 



(•») 



where m* is the typical stellar mass and A includes dimen- 
sionless factors of order unity as well as the Coulomb loga- 
rithm. We use the numerical value of t r given by Binney & 
Tremaine 1987, equation (8-71). Since t r is not very sensitive 
to A we do not attempt to include the dependence of A on 
r. Instead we take 



tr 



1.8 x 10 12 y / 



ln(0.4A) V 100km s" 1 



3 / M Q \ /lO 4 M pc 



,(6) 



where N is the number of stars within the characteristic 
radius 7V 



3.1 The loss cone 

Suppose that stars are removed from the system if they 
come within a distance q of the centre of the cluster. If 
there is a massive black hole at the centre, q will be the 
larger of the tidal radius r T « ii*(M/M*) 1/3 (Hill 1975) or 
the Schwarzschild radius of the black hole. For stars at the 
main sequence turn off, q = tt when M $ 2 x 10 8 A/q. A 
star at radius r will be removed if its angular momentum is 
small enough. Such stars are said to populate a 'loss-cone' 
(Lightman & Shapiro 1977, Young, Shields & Wheeler 1977, 
FR), the angular size of which is given at radius r by 



q 2 q GM 



(7) 



(provided a 1 <C M/q). Even at the edge of the sphere of 
influence, the loss cone is quite small with 8\ c ~ 10~ 7 — 10~ 5 . 

The loss process is usually divided into two regimes ac- 
cording to the angle 8^ through which a star is scattered in 
a dynamical time 



id = r/a. 



(8) 



In the 'diffusive' regime # d < 8\ c , and on timescales much 
longer than t d the loss cone is empty, because a star which 
finds itself in the loss cone is removed within a dynamical 
time. In this case, the capture rate is limited by diffusion 
into the loss cone, and is given, per star, by 
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(9) 



d7V< _ 1 
dN ~~ ln(2/6» lc )t r 

(Lightman & Shapiro 1977), where N(r) is the number of 
stars contained within a radius r. 

In the other regime, 6>d > 0\c> the loss cone is always full. 
Since we assume isotropic velocity dispersions, the fraction 
of stars in the loss cone at any time is just 8\ c 2 . Thus in this 
case the capture rate per star is 



and in case (ii) 



d7V> 
dN 



Pie 



(10) 



The rates (|) and @ are equal at a radius r = r cr it , where 



Gpr A 



M A 



(11) 



where we have written Ai c = ln(2/#i c )/A. The radial depen- 
dence of Aic is weak, so we set it to a constant equal to its 
value at 7V 

The first step in calculating the capture rate is to solve 
equation ( |ll| ) to find r cr i t . Then we can find the total loss 
rate, N, by integrating equations (|^) and ( |lo| ) over the clus- 
ter: 



N = N< + N>, 
with 



N< =4tt 



AT> =4tt 



Ai c G 2 pV dr 
a 3 r 

GMqp dr 



(12) 



(13) 



(14) 



Before moving on, let us examine the various contribu- 
tions to N in some more detail. Let us define /(r) to be the 
r-dependence of the integrand in N < , so 



/(r) = 



G 2 A 3 



Now at small r, a 



-1/2 



(15) 

so /(r) — > as r — > provided 
p < 9/4 (which is true in all the galaxies we consider). The 
lower limit in the first integral in principle should be finite, 
but for p < 9/4 we may set it to zero. At large radius, 

2 2-p 



2 -1/2 

a ~ r ' 



p < 3 
p > 3 



(16) 
(17) 



so f(r) — > as r — » oo provided p > 0. Thus generally /(r) 
will reach a maximum at some radius r max . Quite generally 
we expect that r max will be approximately equal to r a (for 
r a < Tb or < p < 9/4 at large r) or rb (otherwise). We can 
identify two distinct regimes according to whether r cr i t ^ 



(i) If rent S> rmax then N > will be insignificant compared 
to JV<. 

(ii) If r cr it Si r max then A> will be at most comparable 
with 7V< . 



In case |(i)| 



-/V^rnax) 
tr (^max) 



(18) 



N W ATcrit = 



N(r clit ) 



(19) 



tr (^crit ) 

Equations ( |l8| ) and ( p"s| ) constitute a poor man's recipe for 
calculating capture rates, and within a factor of order unity 
are interchangeable with the full integral version in equation 

©■ 

3.2 Red giant capture 

For black hole masses <^ 2 x 10 8 Mq, the Schwarzschild ra- 
dius exceeds the tidal disruption radius for main sequence 
stars. Red giants, as a result of their low densities, have much 
much larger tidal radii, and can therefore be disrupted. Red 
giants may enter the loss cone in the same manner as main 
sequence stars: they diffuse towards an empty loss cone for 
orbits below a critical radius, and above that radius, they 
scatter onto a full loss cone. There is an additional contribu- 
tion to the red giant capture rate. Red giants also grow onto 
the loss cone as they expand. The loss cone for red giants is 
larger than for main sequence stars (q is larger because they 
are less dense) , so even when the main sequence loss cone is 
empty, red giants are being born inside their own loss cone. 
While the total fraction of stars which are red giants at any 
given time may be very small, the number of red giants sus- 
ceptible to capture may be significant. We discuss each of 
these capture channels below. In addition, throughout this 
discussion we use the term red giant loosely to describe stars 
on both the RGB and AGB. 



3.2.1 Radius-time relationship 

In order to calculate the capture rates, it is necessary to 
understand how red giants evolve in radius as a function of 
time because the loss cone, 9\ c oc 7?* (Q). A good estimate of 
the radius evolution can be easily determined as the radius 
and luminosity of a red giant are largely independent of the 
stellar mass, and depend only on the core mass. Functional 
forms for L(m c ) and R(m c ) are given in Joss, Rappaport & 
Lewis (1987) , where m c is the core mass, which are good 
approximations for both the red giant and ascending red 
giant branches. Because the dominant energy source in red 
giants is the p—p chain with ~ 0.7% efficiency and the core 
mass is increased as hydrogen burns to helium, we may write 



L(m c 



O.OO7M 0C 2 f t 



10 



5 V 



-Lr. 



(20) 

i + ro°-v + io°'V iJW ' ^ 

where p — m c /MQ, and the right hand side comes from 
Joss, Rappaport & Lewis (1987). Solving this equation for 
core mass as a function of time, gives the luminosity-time 
dependence. A similar expression to equation ( pl| ) yields the 
radius-time dependence. 

Applying this model between an initial core mass of 
0.17 Mq and a final core mass of 0.8 Mq we find an ini- 
tial radius of 3.09 Rq and a maximum radius, .Rmax, of 
679 Rq after 7 x 10 8 years. This maximum radius is in fact 
larger than the radius predicted by stellar evolution codes, 
because mass loss significantly alters the evolution, espe- 
cially when in the star reaches the thermally pulsing AGB 
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Figure 1. Time as a function red giant radius (solid line). The 
dashed line shows an approximation adopted in eq. 



phase (TPAGB). We therefore limit the red giant evolution 
at 200-Rq which corresponds to the onset of the TPAGB 
(e.g. Bressan etal. 1992). 

For our level of approximation of the red giant capture 
rates, a simple analytical expression for the radius-time de- 
pendence is sufficient. A course approximation (figure [l]) is: 



1 ^0.38 In ( -A 

7xl0 8 years \Rq 



3 + .37. 



(22) 



3.2.2 Red giant birth rate and number density 

In order to calculate the red giant capture rates, we must 
determine their birth rate and their number density. 

The main sequence lifetime of a star of mass m is given 
approximately by 



tn 



10 1 



m. 
~Mr. 







(23) 



and the total number of stars, assuming a Salpeter mass 
function, is 



diV 
dm 



1.6 



M 



© 



(24) 



where the total mass has been normalised to unity in the 
range (0.3, 0.8)Mq. Therefore the birth rate of red giants 
(per star) is 



diVncv 

diV 



AN 
Am 



Am 



At 



0.65 x 10" 



M, 



© 



(25) 



where m is now the main sequence turn off mass (the mass of 
stars which are currently being transformed into red giants). 
The total red giant lifetime £rg ~ 7 x 10 8 y, so the fraction 



of stars that are red giants (with radii between 
200-Rq) is 



/rg 



AN 



tno w 0.05. 



3 and 



(26) 



We therefore expect that the red giant capture rate will be 
lower than the main sequence rate by no more than a fac- 
tor of 20. Another important quantity is the time averaged 
radius, <R> ~ 12Rq. The evolution above WORq makes 
only a small (<10%) contribution to the time averaged ra- 
dius. 



3.2.3 Maximum radius limits 

The red giant capture rates depend on the maximum radius 
reached by the red giant. We consider how this radius may 
be truncated by stellar collisions, and how this radius can 
effectively be truncated when the red giant evolution time 
becomes shorter than the dynamical time. 

First, consider the frequency of collisions between a red 
giant and a main sequence star: 



tcoll 

~t7 



3.71n0.4iV6> 2 



(27) 



l + 26» 

where is the Safronov number, Gm*/2a 2 R (see Binney & 
Tremaine 1987, equation 126). for a ~ lOOkms -1 , 6 varies 
between about 5 and 0.01 as a red giant evolves. The col- 
lisions are important if a red giant typically collides during 
its lifetime: 



At 

tcoll-R 



1, 



(28) 



where -R m i n w 3-Rq. Using the relations from section 3.2.1 
and setting ln(0.4AT) ~ 20, we find that collisions are im- 
portant when 

3 



t r ;S 1.1 x 10' 



^ 250km s" 



years. 



(29) 



For the systems we consider here, the shortest relaxation 
times (at r a ) are many tens of times too large for collisions 
to be important. 

Second, at large system radii, the most extended red gi- 
ants may evolve on a time faster than the dynamical time, so 
that extended red giants have a reduced chance of encoun- 
tering the bl ack h ole. Using the red giant evolution model 



from Section 3.2.1, we calculate the characteristic evolution 
time, R/R as a function of radius (Figure |^). We find that 
over the radii range 3— 200-Rq, the relationship is well ap- 
proximated by 



t R = R/R ~ 6.5 x 10 f 



( R 



\2QQR Q 



(30) 



The maximum radius for a red giant is therefore the smaller 
of 200-Rq and the radius at which tn — id- In practice for 
the Nuker galaxies in the regions of interest tn is always 
much greater than i<j , and we take -R ma x = 200-Rq . 



3.2.4 Capture analogous to main sequence capture 

Main sequence capture occurs by diffusion onto an empty 
loss cone at radii below the critical radius, and by scatter- 
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Figure 2. Characteristic evolution time, R/R as a function radius 
for red giants (solid line). The dashed line shows an simple fit to 
the model. 



ing onto a full loss cone above the critical radius (equa- 
tions ^, In general, the larger the radius of the star, 
the larger the critical radius, and at large radii, the critical 
radius scales as the stellar radius (equations [H], |l^) . For red 
giants we must define three regions because the red giant ra- 
dius and therefore the critical radius are functions of time. 
Below a radius r cr it,i determined by solving equation ( |ll| ) 
with q = 3g ms , the loss cone is depleted on a dynamical time 
for all red giants, so the diffusion rate is appropriate. Beyond 
the second radius r cr i t ,2 determined by solving equation ( |ll| ) 
with q = 200g ms , the loss cone is full, so the full loss cone 
rate is appropriate. Between these two radii, the full rate 
is appropriate until the star reaches a size at which it can 
scatter entirely into or out of the loss cone in one dynami- 
cal time, i.e. the stellar radius Rt for which equation (|ll]) is 
satisfied. The rate of captures of red giants from these three 
contributions can therefore be written: 

diV RG 



dTV 

dTVncw 

dA 



R dt6 lc 
Rq id 



+ 



dt 



Rt 



ln(2/0 lc )t r 



(31) 



The limits in the integrals are understood to mean the time 
at which R is equal to the value specified. The critical red- 
giant radius R t is equal to i? m m for r < r CT it,i and to _R max 
for r > r-crit.2- 



3.2.5 Captures from growth onto the loss cone 

There will be some maximum radius for red giants i? max , 
fixed either by stellar evolution (at ~ 200-Rq), or by other 
means e.g, the rate at which their envelopes are stripped 
by collisions with other stars. As we discuss below, the ap- 
propriate maximum radius is almost always set by stellar 
evolution. This maximum radius defines a loss cone with 
angular size 



r\2 /) 2 Rm&x 



Rr, 



■0 



(32) 



where (as before) 6\ c is the angular size of the main sequence 
loss cone. Assuming that all stars which become red giants 
inside this loss cone will be captured, the rate of red giant 
captures from a given radius (per star) is 



dA RG _ o2 
~dN~ - ° RG 



dN n 



dN 



4606>i c 2 
to 



4.6 x 10 _8 6» lc 2 y _1 , (33) 



where to is the age of the galaxy, and we have chosen R n 
2OOi? . 



3.2.6 Total capture rate 

The rate of red giant capture per star is the sum of equations 
© and df): 



diV 



diV n 



dN 



Rt R dt 6*, c 2 



Rt 



H mm "0 

df 

ln(2/0i c )t r 



Rp, td 



R 



max /) 2 
lc 



Rr. 







(34) 



which must be integrated over the system to obtain the total 
capture rate 



Ar 



47T 



dA RG 2 j 
—pr dr. 



(35) 



For the Nuker galaxies, since r cr it is generally large for main- 
sequence stars it is larger still for red giants, and Rt = Rmin- 
We also find that growth onto the loss cone (equation |33| ) 
dominates the red giant capture rate. 



3.3 Calibration 

FR considered a system in which the stars inside the sphere 
of influence have relaxed to the 'zero-flow' distribution of 
Bahcall & Wolf (1976) with p = 7/4, and outside r a the stars 
are approximately isothermal with a homogeneous core of 
radius ra and density po . One case they concentrated on was 
that in which r cr i t < r a < rrj. CK also calculated the capture 
rate (as well as the full anisotropic distribution function) in 
a system with r cr i t < r a and a relaxed cusp with p ~ 1.8. 
They give a scaling relation (their equation 66) for N in 
these circumstances which is substantially the same as FR 
equation (16a), except for the normalisation. We use this 
normalisation to calibrate our model capture rates. 



4 THE NUKER GALAXIES 

To calculate a capture rate, we need to know the density 
profile p(r) and the black hole mass M. Given the observed 
surface brightness I(r), we can derive p(r) by an Abel inver- 
sion, assuming a constant mass-to- light ratio T. We do not 
assume anything about the density profile in regions which 
are not resolved by the observations, merely extrapolating 
the profiles inwards. This will not affect the capture rates 
provided min[r cr i t , r max ] is not much smaller than the reso- 
lution limit of the observations, which turns out to be the 
case. Most of the data we use are quoted as a Nuker law fit to 
7(r) (Byun et al. 1996) . Thus I(r) is in the form of two power 
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Table 1. Observational data and summary of results. Column (2) taken from Byun etal. 1996; Columns (3-4) from Magorrian etal. 
1998. Column (1), name of galaxy; (2), rb in pc; (3), Nuker law indices (a,/3,7); (4) log 10 M black hole mass; (5) r a /rb; (6) r a in 
arcsec; (7) main-sequence capture time (jV/y); (8) red giant capture time (Nug/y); (9) projected stellar luminosity in central arcsec 
1(1) = L(r < l")/L Edd . 
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Figure 3. Deconvolved V-band surface density of M31 (tri- 
angles) with Nuker law fits to the inner (r < 1.75") and outer 
(r > 3.6") profiles (dotted lines). 



laws, separated by a break radius rb, and p(r) has a simi- 
lar form, with a break radius close to rb. The data we have 
used are reproduced in Table |l|, which also lists some inter- 
esting derived quantities. M32 is listed as NGC 221, M31 is 
NGC 224, and the Milky Way is 'MW. The Milky Way data 
is from the model by Genzel et al. 1996. The Nuker param- 
eters for M32 are taken from Lauer etal. (1992). For M31 



we fit the photometry ourselves using data from Lauer etal. 
(1993) (r < 10") and Kent (1987) (r > 10") and matching 
at r = 10" . The resulting surface density profile is shown in 
Figure [], which also shows the Nuker law fits listed in Table 
0. We fit the nuclear component (r < 1.75") separately, and 
name it NGC 224N in Table |l[ This component is not self 
gravitating and is probably rotationally supported (Lauer 
etal. 1993, Tremaine 1995), and the capture rate is there- 
fore almost certainly much less than in our isotropic model. 
We calculate the latter in order to satisfy ouselves that the 
capture rate is not dominated by the nuclear component, 
which is confirmed by the results. The scale of the flattened 
component in M31 is much smaller than the physical scale 
which dominates stellar capture rates. Similarly, in the other 
galaxies the capture rate is dominated by a physical radius 
which is well resolved, so unresolved substructure should not 
affect our results. 

Throughout this section we assume where necessary 
that the mass of the black hole is proportional to the mass 
of the galaxy (or spheroid) with constant of proportionality 
~ .006 (Magorrian etal. 1998), and the mass to light ratio 
obeys the fundamental plane relation T oc L ' 2 . Thus the 
black hole mass is 



M « 6 x 10 ! 



4.1 



(36) 



Capture rates 

Figure ^ shows the capture rates as a function of the black 
hole mass. In our models all the Nuker galaxies have r cr i t 2> 
rb and all but two have rb <; r a , and hence r max ~ r a . This 
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Figure 4. The capture rates N(y~ ) for main sequence stars 
(squares) and red giants (triangles) as a function of black hole 
mass. The circles are the main sequence capture rates in the case 
where loss cones are all full. The dotted line is the best fit to 
the main sequence rates (with loss cone), and has slope —0.61. 
The dashed line is the upper limit set by Nm* = M/tg with 
to = 10 10 y. 



accounts for why both the main sequence rate and the red 
giant rate are correlated with M. In the red giant case, the 
capture rate is proportional to the mass within r a which is 
of order the black hole mass. In the main sequence case the 
fuelling rate is proportional to this mass divided by tr(r a ). 
The fundamental plane leads to a correlation between t r (r a ) 
and M, and hence to the correlation in Figure [l|. We can 
also calculate 7V max equation (p"8|), and the result is well 
correlated with the value of N given in Table ^] (over a range 
in TV of 3 orders of magnitude TV max is within a factor of two 
or so of TV). We can also use equation jl^ ) to estimate what 
the effect on N would be of errors in the determination of 
the black hole mass. The result for main sequence stars and 
empty loss cones is that errors in M tend to move TV roughly 
along the TV— AT" relation in Figure (dotted line). 



4-1.1 Comparison with supernova rate 

Main sequence capture is faster for smaller black holes, and 
this leads us to expect that the total rate of capture in the 
local universe will be dominated by the galaxies with small- 
est black holes. To see this more clearly, let us compare the 
rate of star capture with the supernova rate in the local 
universe. We assume a galaxy population with luminosity 
function of the form 



dTV 1 / L \ a (-L 



(37) 



(Schechter 1976) where empirically a is typically small and 
negative. We adopt a = -0.07, and L* = 1.8 x 10 10 Lq 
(Estafthiou etal. 1988). The supernova rate in a galaxy is 
roughly proportional to its luminosity L. In the local uni- 
verse it is approximately 



TV SN « 10- 2 iV 1 , (38) 
and the rate of capture of stars from Figure ^ and equation 

(pa is 



(39) 



Together with the fundamental plane relation between T 
and L, Figure § gives s w 0.61 x 1.2 = 0.73. The ratio of the 
star capture rate to the supernova rate is 



/* 



whence 



dTV 
N ^dL dL 
dN 

N S N(L)—dL 



/» ~ 2 x 10 



.4 a+1 / I, 



(—)' 

V Lniin J 



(40) 



(41) 



where Z/ m i n is the smallest luminosity for which we think the 
majority of galaxies have a massive black hole. 

The value of L m i n is not determined by observations, 
since most of the observations are limited to larger systems. 
M32 (L = 3.7 x 10 8 L Q ) is the smallest galaxy in the Nuker 
sample with a confirmed black hole, but it is also an un- 
usually dense system at this luminosity. Dwarf galaxies in 
general are usually thought not to contain black holes. For 
the sake of argument consider L m i n = 10 8 Lq, from which 
equation ^ predicts that for every supernova there would be 
~ 0.01 tidal disruptions. If we take L m i n = lO 7 Z/0 that num- 
ber goes up to ~ 0.1. With searches now detecting many tens 
of supernovae, these numbers suggest that the same surveys 
or similar ones at shorter wavelength should uncover tidal 
disruptions at an interesting rate. 

The rate of red giant captures TVrg is approximately 
proportional to M. From Figure ^] we can read off the con- 
stant of proportionality and use equation (Bfl) to write 



Nr 



10~ 



(42) 



with t ~ 1.2. However, some of the corresponding tidal dis- 
ruption events will have a very long timescale, and hence 
may not be detected as flares. Suppose we can only observe 
flares in cases where the return time for tidal debris t m in is 
less than lOy (see Ulmer 1997 for details). Then the max- 
imum radius of a red giant whose disruption would be ob- 
servable is 

« 1OO7? (^^(JL- y 1/3 . (43) 



1O 6 M 



Since red giant disruption is dominated by stars growing 
onto the loss cone , the rate of disruptions up to i? v is is 
proportional to Rvis, and thus 



iYv 



6 x 10" 



\ 2/3 
tmin \ /_M 

10y I \Lj 



L \ 2t//3 



Integrating over the luminosity function we obtain 

2/3 

-ii- ', 

/.vis « 6 x 10 



'mm 

lb? 



a + 1 



2£ + 3a ' 



(44) 



(45) 



where /* v is is the ratio of the rate of visible red giant dis- 
ruptions to the supernova rate, so the red giant capture rate 
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is probably much less than the main sequence capture rate 
in the local Universe. Note however that it is not possible 
to rule out intermittent behaviour in the accretion of tidal 
debris (e.g. Lee, Kang, & Ryu 1996) and so equation ( plo] ) 
represents a lower limit to the rate of detectable red giant 
disruptions. 

4-1-2 Upper limits to capture rate 

An upper limit to the capture rate is often given as the full- 
loss cone rate in the region r < r cr i t . If non-axisymmetric 
processes can keep the loss-cone full the capture rate can 
be approximately derived from equation (fid). The rate of 
capture of stars from the full loss cone would be given by 
N > but integrated all the way down to r = r a (assuming 
that the black hole itself imposes enough symmetry to keep 
the loss cone starved at r < r a ). The dominant radial scale 
is thus that at which p/a is largest. Since a ~ r -1 / 2 for 
small r, we see that if p ~ r~ p then N is dominated by 
the contribution from r a if p > 1/2. This is the case in our 
models in all but three cases (and a more realistic model of 
the bound stars inside r a would probably also havep > 1/2). 
The full loss cone rates for the Nuker galaxies are of order 
10~ 3 y~ 1 , roughly independent of M (circles in Figure |^). 
Thus the capture rate would be roughly 

A/to* = M/to ~4x 10~ 4 (J^) M Q y-\ (46) 

with s « 0.2, which is larger than the supernova rate for L 
less than a few percent of L* . Substituting equation (^) in 
equation (|4o|), we find 

/.« 4 x 10- 2 SL±1 ( hsz. Y~ a (47) 
s — a V L* / 

where L max is the largest luminosity galaxy which can tidally 
disrupt a star (as opposed to swallowing it whole). Tak- 
ing the maximum black hole mass for tidal disruption to 
be 10 8 Mq we obtain /* ~ 0.1 . This means that tidal dis- 
ruptions would be about ten times less likely than super- 
novae if the loss cone was full. The tidal disruptions in 
this case would also be coming mainly from galaxies with 
L ~ Z/nmx- The red giant capture rate would be simply 
/RG-Rm ax /i?0 ~ 10 times the main sequence rate (with 

Equations (^) and ( fi"l| ) together provide a constraint 
on the mass supply for the black holes. If only a small num- 
ber of tidal disruptions were observed, and they came from 
the faintest galaxies, then the dominant fuelling would be 
via empty loss cones. 

Another upper limit to the capture rate is given by as- 
suming that the galaxies we see today are not in a special 
state, so their capture rates must be smaller than M/to (the 
dashed line in Figure ^]). In fact this is a natural value for 
the capture rate in the case that the black holes grew by 
eating stars from the cluster around them. For the majority 
of the Nuker galaxies M/to is actually larger than the full 
loss cone rate. Thus there has not been enough time for the 
black holes to grow by eating stars from the nuclear star 
cluster. A more efficient mechanism must be responsible for 
their having grown to the present masses. Merrit & Quinlan 
(1997) have argued that the universal x ~ 0.006 is a result 
of the black hole eating stars until the intrinsic triaxiality of 
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Figure 5. The radius of the sphere of influence r a in arcsec as 
a function of black hole mass. 

the star cluster is reduced by the presence of the black hole 
(Gerhard & Binney 1985). At least for the larger galaxies, 
triaxiality could never supply enough stars to provide all the 
mass in the black holes. 



4.2 Keplerian regime 

The criticism that r a is often comparable with the resolution 
limit of the observations in question (Rix 1993) is no longer 
valid. This has been pointed out by Kormendy (1994), and 
we support his view with our carefully calculated values of r a 
in the sense that the measured values do not cluster around 
any particular value (Table [j] and Figure |B|). On the other 
hand, some of the smallest black holes have r a comparable 
to or smaller than 0.1". 

The black hole masses measured by Magorrian etal. 
(1998) required detailed modelling of the kinematics of the 
surrounding galactic nuclei. In principle one should see at 
small radius the Keplerian regime in which a 2 ~ M/r. This 
would be a very clear signature of a black hole, giving di- 
rectly the mass to within a factor of order unity. The Kep- 
lerian regime should occur within some radius of order r a . 
We find that the logarithmic slope of a 2 reaches —0.95 at a 
radius tk which is generally much less than r a , with median 
value ~ 0.2r a (this is of course consistent with the fact so 
much modelling had to be done to measure the black hole 
masses) . This situation is also seen in the centre of the Milky 
Way, where the Keplerian regime is now resolved (Genzel 
etal. 1996) at r K ~ r a /10. 

One might expect that the galaxies with the largest val- 
ues of 7~k would have the most secure black hole detections, 
but this appears not to be the case. The three galaxies where 
we find that r K > 0.5" are NGC4486B, NGC7768 and M31, 
and these do not appear to have significantly better black 
hole masses than any of the others. 
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5 DISCUSSION 

Note that the scaling relations of FR and CK giving ./V 
in terms M, po and ro apply only in very specific circum- 
stances. The result of CK only applies to the case r cr i t < r a . 
If fcrit < r a and p ~ 1.8 for r < r a then CK equation 
(66) can be used to calculate N, replacing po — > p(r a ) and 
ro — ► r a , with the result roughly independent of what hap- 
pens at r > r a . FR are careful to distinguish between the 
cases r cr it ^ r a , and in these cases they use p — 7/4 or 
p — accordingly. It is straightforward to generalise their 
expressions using equations ( |l9| ) and ([bsl). 

Rauch & Tremaine (1997) have discussed the enhance- 
ment of tidal disruption rates by a process they call 'resonant 
relaxation'. This effect can increase tidal disruption rates of 
stars which are bound to the black hole (i.e. at r < r a ) by 
factors of order unity. Detailed calculations by Rauch & In- 
galls (1997) also show that the effect is also quencehed by 
relativistic precession for black hole masses M <; 10 8 Mq. In 
the models we constructed the fuelling rate was never dom- 
inated by strongly bound stars, so in common with Rauch 

6 Ingalls (1997), we find that resonant relaxation has little 
effect on tidal disruption rates. 

There have been some observations of active galaxies 
which develop a broad-lined feature in their spectrum over a 
period of the order of months (Storchi-Bergmann et al. 1995, 
Eracleous etal. 1995). Eracleous etal. (1995) constructed a 
phenomonological model of such an event which involved 
an elliptical accretion disc, and suggested that such a disk 
could arise from the tidal disruption of a star (see also Syer 
& Clarke 1992). The size of the disc they required to fit the 
observations was too large to correspond to a main sequence 
disruption. It could, however, arise f rom the disruption of a 
red giant. Our discussion in Section 4.1 concluded that red 
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giant disruptions with such short timescales should be very 
rare, so it is perhaps not surprising that such events are not 
encountered frequently. 
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